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ABSTRACT 

Little is known about the incidence of magnetic fields among the coolest white dwarfs. 
Their spectra usually do not exhibit any absorption lines as the bound-bound opacities 
of hydrogen and helium are vanishingly small. Probing these stars for the presence of 
magnetic fields is therefore extremely challenging. However, external pollution of a 
cool white dwarf by, e.g., planetary debris, leads to the appearance of metal lines 
in its spectral energy distribution. These lines provide a unique tool to identify and 
measure magnetism in the coolest and oldest white dwarfs in the Galaxy. 

We report the identification of 7 strongly metal polluted, cool ( T e g < 8000 K) 
white dwarfs with magnetic field strengths ranging from 1.9 to 9.6 MG. An analysis of 
our larger magnitude-limited sample of cool DZ yields a lower limit on the magnetic 
incidence of 13 ± 4percent, noticeably much higher than among hot DA white dwarfs. 

Key words: stars: white dwarfs - stars: magnetic field - stars: planetary systems - 
stars: evolution 


1 INTRODUCTION 


White dwarfs (WDs) have been known to harbour mag¬ 
netic fields si nce the detection o f circularly polarised light 
from GJ 742 jKemp et al.lll970l ). In the following decades 
a plethora of magnetic WDs (MWDs) have been iden¬ 
tified either from Zeeman splitting of abs orption lines 
in th eir spectra or by spectropolarimetry jKawka et al.1 
12001 and references therein). A wide variety is seen in 
temperature, atmospheric composition, and field strength. 
The advent of large scale spectroscopic surveys, in par¬ 
ticular the Sloan Digital Sky Survey (SDSS), has in 
the last decade increased the number of known MWDs 


2003i: Vanlandingham et al. 

20051; Kleinman et al.1 

Kepler et al.ll20li 120151). 



Despite the ever growing list of these previously rare 
objects, two questions continue to remain without a definite 
answer: What is the origin of these magnetic fields? And 
what is the fraction of WDs that are magnetic, and how 
does this vary with cooling age/temperature? 

Two distinct models have been proposed to explain 
the emergence of fields > 1 MG in isolated WDs. In the 
fossil field hypothesis, the magnetic fields of the chemi¬ 
cally peculiar Ap/Bp stars are thought to be amplified 
due to flux conservation during post-main sequence evo- 
lution resulting in WDs with fields in the MG regime 
JWoltjei 196^_ Angel fe Landstreetlll97d : lAngel et al.lll98ll : 
IWickramasinghe fe Ferrariol 200d) . A more recent hypoth¬ 


esis dTout et al.l 120081 ) considers a binary origin, where a 
system undergoing a common envelope leads to magnetic 
dynamo generation. 

The incidence of magnetism in WDs remains poorly es¬ 
timated due to selection effects. Independent studies are dif¬ 
ficult to reconcile with one another as each suffers from its 
own set of biases. This problem becomes significantly more 
pronounced when focusing on subsets of the total WD pop¬ 
ulation where small number statistics dominate. Recent vol¬ 
ume limited samples of nearby WDs present the most unbi¬ 
ased estimates of the magnetic incidence when considering 
all WD sub types, and suggest incidences of 21 ± 8 percent 
for WDs wit hin 13 pc of the Sun , and 13±4 percent for those 
within 20 pc dKawka et al.ll2007l ). However these MWDs are 
dominated by fields lower than 100 kG and strongly mag¬ 
netic objects with fields above 10 MG. On l y 1 ou t of the 15 
MWDs in the compilation of lKawka et al.l (120071 ) has a field 
strength between 1 and 10 MG (the range th at we discuss in 
this work). More recently, ISion et all J2014I ) have presented 
a volume limited WD sample within 25 pc from the Sun. 
They find a magnetic incidence of 8 percent when consid¬ 
ering magnetic fields above 2 MG only. Other studies have 
investigated the magnetic incidence with much larger, but 
magni tude-limited samples. For instance iKIeinman et al.l 
(l2013l) identified over 12000 DAH from SDSS data release 7 


1 WDs showing only hydrogen/helium lines in their spectra are 
classified DA/DB, with only metal lines as DZ, and without any 
spectral lines as DC. Magnetic DA WDs where magnetism is de- 
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(DR7 ) spectra, of which over 500 are suggested to be mag¬ 
netic dKepler et al.ll2013l) . leading to a much lower incidence 
of 4percent. However, because this sample is magnitude- 
limited, it is intrinsically biased. Most degenerates in the 
local sample have tempera tures below 10000 K, whereas 
84 percent of the WDs from iKepler et al.l (120131 4 are hotter 
than this. The discrepant numbers between the local sample 
of cool/old WDs and hotter/younger WDs, have been the 
basis for some auth ors to claim an age-dependency of the 
magn etic incidence dFabrika fe Valvavinlll999l : lLiebert et al.l 
1 2003li . 

An alysing t h e sma ll sample of WDs with accurate par¬ 
allaxes, iLiebertl d 19881 4 noted that magnetic WDs appear 
to be under-luminous for their colour, suggesting they have 
smaller radii, and hence hi g her m asses, than non-magnetic 
WDs. Later, iLiebert et aid d2003d ) derived a mean mass of 
0.93 M 0 for eight MWDs from the Palomar Green (PG) 
survey, based on model atmosphere analyses, compared to 
~0.6 M@ for non-magnetic ones. While there is hence inde¬ 
pendent evidence for higher-than-average masses for MWDs, 
caveats to bear in mind are that there are still few MWDs 
with precise parallaxes, and even for those systematic uncer¬ 
tainties in the analysi s of their spectra l imits the accuracy 
of the desired masses dKiilebi et al.ll201ol ). 

A common theme among all the above investigations is 
that the true magnetic incidence is expected be higher, as 
the various biases (e.g. signal-to-noise, magnetic broaden¬ 
ing) tend to work against the identification of MWDs. 

The very coolest WDs (T e ff< 8000 K) do not show 
optical lines of hydrogen or helium since these elements 
are in their ground states in the low temperature atmo¬ 
spheres. This transition to featureless (DC) spectra occurs 
at around 11000 K for WDs with helium dominated atmo¬ 
spheres, an d 6000 K for WDs with hydrogen dominate d at¬ 
mospheres (IVauclair et al.lIl98ll ; iBergeron et al.l l200ll) . Be¬ 
cause of this absence of absorption lines, it is not possible 
to identify magnetism in these stars via Zeeman splitting. 
A handful of DCs h ave been foun d to be magnetic through 
spectropolarimetry (|Putnevlll997l ). However, this method is 
expensive and is unsuitable for most known cool WDs be¬ 
cause of their faintness. 


2 METAL POLLUTED WHITE DWARFS 

The last two decades have seen great interest devoted to 
the study of WDs with atmospheres contaminated by met¬ 
als. The now widely accepted scenario is that after post- 
main sequence evolution to the WD stage, an accompa- 
nying planetary syste m will undergo dy nam ical instability 
dDebes fe Sigurdssonl 120021 ; IVeras et al.l 120131 ). Small rocky 
bodies such as asteroids or minor-planets may then have 
their orbits perturbed by a larger planetary object and 
soon find themselves venturing into the Roche-radius of the 
WD, resulting in their tidal disruption dDebes et al.lI20l3 : 
IVeras et al.l 120141 ) and subsequent formation of a circum- 
stell ar debris disk. This mate rial is then accreted onto the 
star djurall2003l : lRafikovll201ll ) producing metal lines in the 
spectra of these WDs. 


tected via Ze eman splitting ar e known as DAH, and via polarime- 
try as DAP dSion et al.lll983l ). 


Co ol WDs displaying only metal lines are classified 
as DZ. iFarihi et aid J2OI0I ) showed that DZ and DC white 
dwarfs share the same velocity, spatial and temperature dis¬ 
tributions and should therefore belong to the same stellar 
population. Assuming the existence of magnetism among DZ 
stars is uncorrelated with the presence of metals, the detec¬ 
tion of split metal lines in the spectra of DZ white dwarfs 
becomes a powerful tool for determining the magnetic field 
strength distribtion among the very coolest and oldest of 
WDs. However, prior to t his wor k, onl y 3 m agnetic. DZ WDs 
were known: LHS2534 dReid et ah 200ll). W D0155+003 
dSchmidt et al.ll2003l ). and G165-7 ( Dufour et alfeOOtil) . with 
respective surface averaged field strengths, Bs, of 1.9, 3.5, 
and 0.6 MG. These all have SDSS spectra and so we re¬ 
fer to these using the SDSShhmmiddmm naming format 
(SDSS1214—0234, SDSS0157+0033, and SDSS1330+3029 
respectively) for consi stency with our new identifications 
and the name s used inlK oester et aid d201ll) . 

Recently, iKoester et al. ( 20 111 ) identified a sample of 
26 cool (T e fi< 9000 K) DZ with strong photospheric metal 
pollution, filling previously empty parameter space at low 
T e ff and high atmos pheric Ca abu nd anc e compared with the 
preceding work by iDufour et all d2007t) . The rocky nature 
of the accreted material is evident from the variety of de¬ 
tected metals which includes Ca, M g, Fe, Na, Cr, an d Ti, 
and the low abundance of H. The IKoester et ahl d201ll ) sam¬ 
ple also includes re-identifications of the known magnetics 
SDSS0157+0033 and SDSS1330+3029 (see above). 


3 RESULTS 
3.1 Identification 

In an extension of the work by IKoester et all (l201ll ) using 
the more recent DR10 spectroscopic data, we have iden¬ 
tified a sample of 79 strongly polluted cool DZ, with the 
serendipitous discovery that 10 are magnetic. Full details of 
this sample will be presented elsewhere (Hollands et al. in 
prep); here we only provide a brief summary of the selection 
procedure. _ 

Firstly, we queried SDSS Casjobs (iLi fc Thakail 120081) 
for objects with colours below the main sequence in the 
(u — g) vs. (g — r ) two-colour diagram (Figure [T]l. This re¬ 
gion of colour-space is highly abundant in quasars of redshift 
2-5. To remove these, we employed two independent cuts 
in proper-motion and redshift: a minimum proper-motion 
threshold of 3-cr above zero and SDSS spectroscopic redshift 
below 0.01. Any spectrum with non-zero redshift zwarning 
(0 indicates no warning) was also kept. 

Requiring the objects to pass at least one of the proper- 
motion or redshift selection criteria was necessary as neither 
could be completely relied upon in isolation. On the one 
hand, proper motions may be unavailable if the object is too 
faint to be detected in the USNO-B images used to detect 
positional variations. Additionally, because detectable cool 
DZ are by definition relatively nearby (as they are intrinsi¬ 
cally faint), some will have sufficiently high proper-motions 
to evade a match between SDSS and USNO-B observations. 
On the other hand, the unusual spectra of DZ WDs with 
their deep, broad metal lines can trick the SDSS redshift 
pipeline into classifying them as quasars, reporting values of 
Z > 1, with ZWARNING = 0. 
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Figure 2. Optical SDSS spectra for all 10 magnetic DZ, ordered from top to bottom by increasing magnetic field strength. Data 
are smoothed by a 7-point boxcar to reduce the noise level and enhance spectral features, with the exception of the bright object 
SDSS1330+3029. Insets are plotted where splitting is not obvious in the full spectrum (same smoothing applied). 
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Figure 1. The colour-cut used to identify cool DZ WDs extends 
below the main sequence (dashed line). Previously known mag¬ 
netic DZ are indicated by filled triangles, and new discoveries by 
filled circles. Crosses indicate the other 69 DZ in our sample of 79 
which do not show Zeeman splitting of spectral lines. The main 
sequence is indicated by the grey-scale. 


In addition to quasars, low mass main sequence stars 
strayed into our colour cut and were removed via template 
fitting for K and M spectral types. Finally, the DZ WDs were 
identified by visual inspection of the remainin g spectra. This 
metho d recovered all 17 DZ WDs presented bv lKoester et al.1 
(1201 ll) that fell into our colour-cut, while also finding 62 
additional objects, bringing our total colour-selected sample 
to 79. Magnetic DZ were identified from splitting of the Mg I 
blend at 5172 A and the Nal doublet at 5893 A. 

In total we found our sample of 79 DZ contains 10 
MWDs, which are listed in Table [T] with their SDSS 
optical spectra shown in Figure [2] Among these 10 
objects, 3 are the previously known magnetic DZ as 
mentioned in M (SDSS0157+0033, SDSS1214-0234, and 
SDSS1330+3029). One s ystem, SD SS11 5 2+16 05, was first 
identified as a DZ by iKoester et al.l d201 il l . who also 
pointed out the possible magnetic nature of the star. How¬ 
ever the signal-to-noise (S/N) of the single SDSS spec¬ 
trum available was too low for a firm conclusion. The 
newer SDSS spectrum, presented here, shows unambigu¬ 
ous Zeeman splitting. The remaining 6 objects are all 
newly identified WDs. We note tha t the recently pub¬ 
lished WD list of iKepler et all (l2015l i independently iden¬ 
tifies SDSS0735+2057, SDSS0832+4109, SDSS1003-0031, 
and SDSS1536+4205, however only SDSS1536+4205 is re¬ 
ported to be magnetic. 


3.2 Average magnetic field measurement and 
effective temperatures 

In the spectra of all 10 objects (with the exception of 
SDSS1330+3029 where Zeeman splitting is only seen in the 


Ca triplet), we are able to identify split lines of either Mg I, 
Na I, or both. For field strengths of > 2 MG, spin and orbit 
angular momenta decouple, and so the Paschen-Back ap¬ 
proximation is appropriate. Therefore we treat the splittings 
as triplets. 

We used a 7-parameter fit to the observed line profiles 
to measure the field strength. The continuum flux in the 
vicinity of the triplet was modelled as a quadratic in F v . 
A linear approximation would not suffice, particularly for 
the wings of the broad Mg feature (see Figure [2j. We then 
modelled the triplet as the sum of three Gaussian profiles 
with equal width, depth (in continuum normalised flux) and 
separation in wavenumber, 1/A. The wavenumber of the n- 
component of the triplet was also included as a free param¬ 
eter to account for small shifts. In all cases, we found the n- 
components are blue shifted from their rest wavelengths, and 
generally increasing with field strength, suggesting this is 
predominantly caused by the quadratic Zeeman effect, with 
only minor contributions from gravitational and Doppler 
shifts. The maximum blueshift of the 7r-component is found 
to be 5 A for the Na triplet of SDSS1536+4205. The small 
(few percent) measurement error this may have on our field 
measurements does not affect our discussion on magnetic 
incidence. 

We used least-squares minimisation via the Levenberg- 
Marquardt algorithm to optimise these parameters. Where 
possible we fitted both the Mg and Na lines, however, this 
could not always be achieved for a variety of reasons: One of 
the lines may be significantly less deep than the other; the 
Mg line in some cases is very broad and asymmetric, such 
that the 3 components cannot be distinguished; or poor sub¬ 
traction of sky emission distorts the flux near the n compo¬ 
nent of the Na triplet, making a fit to this line less reliable 
than for the Mg triplet. 

The average surface magnetic field strength, Bs, was 
subsequently calculated from 


B s / MG 


A(l/A) 
46.686 ’ 


(1) 


where A(l/A) is the inverse wavelength se paration in 
cm" 1 between the components of a triplet dReid et all 
l200lf) . As an example, the fit to the most magnetic object 
SDSS 1536+4205 is shown in Figure[3]with a measured field 
strength of Bs = 9.59 ± 0.04 MG. 

Additionally, we estimated T e ff, by fitting each SDSS 
spectrum to a grid of DZ model spectrtQ The grid 
was calculated for temperatures from 4400 to 14000 K in 
steps of 200 K, and Ca abundances from log 10 [Ca/He] = 
— 10.50 to —7.00 in steps of 0.25 dex. Because photospheric 
metal lines arise from the accretion of pl anetary debris 
(IZuckerman et akll2007l : rGansicke et alj|2012h . all other met¬ 
als were held at a fixed ratio to Ca, assuming bulk Earth 
abundances. The models were calculated for fixed surface 
gravity of logg = 8.0. 

The SDSS spectra were fit to each of the model spectra 
to obtain a \ 2 , at each grid point. The x 2 {T e a, log 10 [Ca/He]) 
plane was then interpolated with bicubic splines to more 
accurately locate the best fit in the parameter space. The 


2 As our non-magetic models ignore the potential additional line 
blanketing from the Zeeman-split lines, the true effective temper¬ 
atures may differ from the values quoted here 
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Table 1. Magnetic DZ with field strengths, temperatures and SDSS PSF photometry. Previously k nown magnetic D Z ar e indicated by 
their starred coordinates. The field strengths for SDSS1214—0234 and SDSS1330+3029 are taken from lReid et alJ (j 200 il l and lDufour et all 
(POOfl ) respectively. The quoted By uncertainties are obtained from the formal errors on our fit parameters. Where a good fit is obtained 
for both Mg and Na lines, their weighted mean is used. 


SDSS 

Bs (MG) 

T eff (K) 

u (mag) 

9 (mag) 

r (mag) 

i (mag) 

z (mag) 

J133059.26+302953.2* 

0.65 

6000 

18.28 + 0.02 

16.32 + 0.02 

15.89 ± 0.01 

15.91 ± 0.02 

16.09 ± 0.03 

J090222.98+362539.6 

1.92 + 0.05 

6300 

22.58 + 0.32 

20.75 + 0.03 

20.46 ± 0.04 

20.36 ± 0.06 

20.50 + 0.22 

J121456.39—023402.7* 

1.92 

5200 

20.87 + 0.06 

18.32 + 0.02 

17.75 ± 0.01 

17.56 + 0.01 

17.51 ± 0.02 

J083200.38+410937.9 

2.35 + 0.11 

5900 

23.62 + 0.79 

21.32 + 0.05 

20.63 ± 0.03 

20.70 + 0.05 

20.75 ± 0.14 

J115224.51+160546.7 

2.72 + 0.04 

6500 

21.73 + 0.12 

20.18 + 0.03 

19.95 ± 0.02 

20.02 ± 0.03 

20.08 + 0.09 

J015748.14+003315.0* 

3.49 + 0.05 

5700 

21.30 + 0.07 

19.59 + 0.02 

19.19 + 0.02 

19.21 ± 0.02 

19.36 + 0.05 

J100346.66—003123.1 

4.37 + 0.05 

6300 

22.91 + 0.33 

20.61 + 0.03 

20.06 ± 0.02 

20.00 ± 0.03 

20.20 + 0.11 

J073549.19+205720.9 

6.12 + 0.06 

6000 

23.08 ± 0.34 

20.53 + 0.02 

19.98 + 0.02 

19.92 ± 0.02 

20.09 ± 0.09 

J232538.93+044813.1 

6.56 + 0.09 

7200 

21.44 + 0.15 

19.88 + 0.02 

19.58 + 0.03 

19.68 + 0.04 

19.69 + 0.11 

J153642.53+420519.2 

9.59 + 0.04 

5500 

23.32 + 0.65 

20.84 + 0.04 

20.30 + 0.03 

20.17 + 0.04 

20.33 + 0.16 



Figure 3. Fits to the Mg and Na splittings for SDSS1536+4205. 
Lines are fit with Gaussians with equal 1/A separations from the 
central rr components. 


corresponding effective temperatures are quoted in Table [T] 
to the nearest 100 K. 


3.3 Field measurement for an offset dipole 

For a MWD with the simplest possible field structure, a 
centred magnetic dipole, the magnetic field across the star 
varies with magnetic latitude, resulting in a field twice 
as strong at the pole s compared to the magnetic equator 
(lAchilleos et al.lll992l ). A spectrum taken of a MWD inte¬ 
grates over its entire visible surface, and therefore over the 
range in field strengths, resulting in magnetic broadening of 
the a components in a given Zeeman triplet. This effect is 
generally observed in magnetic WDs and can be used to con- 
strai n the inclination to t he magnetic axis to some degree 
(e.g. iBergeron et al.|[l992l ). 

A brief inspection of Figure [2] shows that some of the 
MWDs in our sample have very narrow Zeeman triplets, 
where all three components have similar depths. This is seen 
in SDSS0157+0033, SDSS1003-0031, SDSS2325+0448, and 
in particular SDSS1536+4205. A centred dipole field ob¬ 


served from any inclination to the magnetic axis can not 
reproduce these line profiles. 

This suggests that these WDs may have more complex 
field topologies. We show that a dipole offset from the star’s 
centre can reproduce the observed Zeeman line profiles. In 
principle this offset, a, can be i n any direction relative to the 
unsh ifted magnetic field axis dAchilleos fc Wickramasinghel 
Il98dh . however, here we consider only displacement along 
the magnetic dipole axis (whic h we define to be in t he z- 
direction), i.e. a x = a y = 0 as in lAchilleos et abl (Il99 2tl. 

For an arbitrary point on the surface of the WD with 
coordinates (x, y, z) in unit s of 7?wd, the stren gth of the 
field, B(x,y,z), is given by dAchilleos et al.lll992l ) 

B(x, y , z) = B d [r 2 + 3 (z - a z ) 2 ] 1/2 /2r 4 , (2) 

where a z is the dipole offset, B d is the dipolar field strength, 
and 

r 2 = x 2 + y 2 + {z - a z ) 2 . (3) 

These equations can be used to compute synthetic Zeeman 
line profiles given for a given B d , a z , and line-of-sight incli¬ 
nation to the magnetic axis, i. In our model we firstly gen¬ 
erate a set of 10000 points randomly distributed over the 
projected surface of the WD. At each point the magnetic 
field strength is evaluated using equation ©, accounting 
for the inclination of the magnetic axis, and the correspond¬ 
ing Zeeman line profiles are generated as a sum of 3 Gaus¬ 
sians separated according to equation ©. These profiles are 
then coadded with weights proportional to the limb darken¬ 
ing corresponding to the location on the star. We use limb 
darkening coefficients app r opriat e for a 6000 K, log <7 = 8 
WD from iGianninas et al.l (|2013l ). adopting the logarithmic 
limb darkening law described therein. 

We fitted the above model to the Zeeman lines in 
SDSS 1536+4205 using the affine invariant MCMC sampler, 
emcee. Replacing Bs with B d and including the inclination 
and dipole offset, increased the number of free parameters 
in the fit (compared with 1 13.21 to 9. Uninformed priors were 
used for variables with a physically constrained range. E.g. 


3 B d is defined in such a way that the magnetic field has this 
strength at z = a z ± 1 WD radii along the magnetic axis. 
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Wavelength [A] 


Figure 4. Similar to Figure [3] but fitted with an offset dipole 
model. The adopted values of B d , i and a z are given in Table [2] 
i.e. the two triplets are shown with their individual fits. The a 
components are now slightly less deep than the respective tc com¬ 
ponents, as seen in the data. 


— 1 < a z < +1 forces solutions with the dipole centre con¬ 
fined within the stellar surface. For the line-of-sight inclina¬ 
tion to the magnetic axis, the prior distribution P(i) oc sin i 
was employed. 

Since SDSS1536+4205 shows distinctly split lines of 
both Mg and Na, both were fitted independently. The values 
for Bd, i, and a z are shown in Table [2] with the resulting 
best fits to the spectra shown in Figure]!] While the line pro¬ 
files appear similar to the simple model shown in Figure [3] it 
should be recalled that we have now fitted a physical model 
capable of reproducing the observed narrow Zeeman lines 
rather than the assumption of unbroadened lines used in 

The fit to the Mg triplet has a slightly worse reduced 
X 2 than to the Na triplet. This is the result of the Mgl 
line having an intrinsically asymme tric profile due to quasi¬ 
static broadening o f this transition (IWehrse fc Lieber3ll980l : 
iKoester et al.luOllI '). Therefore assuming a Gaussian profile 
limits the quality of the fit. Nevertheless, the resulting pa¬ 
rameters from the Mg and Na fits are in agreement within 
their (similarly large) uncertainties. 

While the inclination uncertainties permit a wide range 
of values within the allowed parameter space (0-90°), the 
results for Bd and a z strongly suggest an offset dipole. How¬ 
ever, it should be noted that Bd and a z are highly anti¬ 
correlated, with a correlation coefficient of —0.99 for both 
fits. The values of a z we find a re well within the range of 
those found for SDSS DA WDs (iKfilebi et al.ll2009h . 

The fit values of Bd, i, and a z have the following phys¬ 
ical interpretation: An offset dipole leads to a strong field 
emerging at one of the poles of SDSS1536+4205, (fa 50 MG 
according to equation ©), with the opposite hemisphere 
exhibiting a very uniform field strength of 9.6 MG as in Ta¬ 
ble [T] The value of inclination and sign of a z imply that 
most of the WD surface with high fields is obscured when 
viewed from the Earth, with only a small amount of this re¬ 
gion entering the limb of the star. Hence the a components 


Table 2. Results from our MCMC fits to SDSS1536+4205 as¬ 
suming an offset dipole field structure. Quoted values and uncer¬ 
tainties errors are, respectively, the 50tli, 15.9/84.1th percentiles 
of the posterior probability distributions. Reduced \ 2 values are 
calculated using the median for each parameter. 


Parameter 

Mg 

Na 

B d (MG) 

18.6ti| 

20.0+j 

inclination (°) 

CO 

M■ 

1 + 

cn ^ 

41+3-6 

a z (%%d) 

-23.9)+f 

-28.8j)4 3 

reduced x 2 

1.25 

0.97 


of the Zeeman lines are broadened only slightly with their 
depths reduced by a few percent. 


4 NOTES ON INDIVIDUAL OBJECTS 

SDSS0157+0033. The value of Bs given bv ISchmidt et all 
(120031') is 3.7 MG, whereas we obtain a value of 3.49 ± 
0.05 M G using the same SDSS spectrum. ISchmidt et al.l 
(|2003t) arrive at this value from measuring both the Mg and 
Na lines. For the Na triplet, only two of the Zeeman com¬ 
ponents can be identified in the noisy spectrum, and their 
centres are difficult to locate. For the Mg triplet, where all 
three Zee man component s are well resolved, we have no¬ 
ticed that ISchmidt et al.l (120031 ') report a position for the 
cr~ line (5128 A) which is about 7 A bluer than we mea¬ 
sure leading to an overestimate of Bs We therefore suggest 
3.49 ± 0.05 MG as a revised value for this star’s surface av¬ 
eraged magnetic field strength, which we quote in Table [I] 
SDSS0735+2057 exhibits an extremely broad Mg feature, 
and close examination reveals Zeeman split lines at its base. 
None of the other 78 DZ in our sample shows an Mg line like 
this. It seems likely that the unusual Mg profile is a direct re¬ 
sult of the magnetic field partially splitting an already broad 
feature. SDSS2325+0448 has a similar surface-averaged field 
strength and so one might suspect a similarly broadened Mg 
triplet, which is not observed. The stark difference between 
their line profiles can be explained by higher metal abun¬ 
dances and lower T e ff for SDSS0735+2057. This increased 
opacity that the magnetic induces in the Mg feature will 
need to be considered when calculating chemical abundances 
(Hollands et al. in prep). The formal errors on Bs (Table]!]) 
are surprisingly small, at the level of ~1 percent, however 
the values obtained from the Mg and Na lines for this star 
agree to within 1.3cr (6.24 ± 0.11 MG and 6.08 ± 0.07 MG 
respectively). We adopt the weighted mean of these as the 
measured value in Table [T] 

SDSS1536+4205 has the strongest magnetic field of any 
known DZ, with Bs = 9.59 ± 0.04 MG. We also obtain very 
consistent field measurements between the Mg and Na lines, 
9.57 ± 0.08 MG and 9.60 ± 0.04 MG respectively, indicating 
the uncertainties are not underestimated. Again the final 
value presented in Table [l] is the weighted average of these 
two independent measurements. The Zeeman split lines of 
Mg and Na for this WD are well resolved from the noise, 
yet do not show any significant magnetic broadening, as 
discussed in detail in £ 13.31 and H5.5I We also note that this 
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object may show an emission line at approximately 8400 A. 
This peculiar feature does not coincide with any sky line, 
and peaks at 4cr above the continuum. Additionally, the fea¬ 
ture is visible in multiple SDSS sub-spectra. Checking spec¬ 
tra from adjacent fibers observed on the same night and 
plate does not indicate flux contamination from other fibers, 
and neither does inspection of the SDSS images reveal any 
nearby bright stars. If this emission is real, the most plau¬ 
sible identification is Ol. W e note that a few cool DQ WDs 
show oxygen emission lines dProvencal et al.ll2005l ). the ori¬ 
gin of which is still uncertain. A higher S/N spectrum will 
be required to determine whether this feature in SDSS is 
genuine. 


5 DISCUSSION 


5.1 Magnetic incidence 

We have amassed a sample of 79 DZ WDs with 
colours that place them below the main sequence in the 
(u — g ) vs. (g — r) plane. We identify 10 of these to be 
magnetic (Table [I]), leading to an observed incidence of 
13 ± 4 percent purely for this sample. Our selection proce¬ 
dure only uses colour, redshift and proper-motions. Since 
colours will not be significantly altered by the presence of 
a magnetic fiekfl we believe this selection procedure to be 
unbiased towards (or against) selection of magnetic objects. 
This is not to say this sample is free from biases. On the 
contrary, various selection effects, which we discuss in the 
following sections, suggest that our measured magnetic in¬ 
cidence is only a lower limit. 

We report DZ to be magnetic only where we confidently 
detect Zeeman splitting. The major limitation in detecting 
splitting of spectral lines are S/N and spectral resolution. 

The vast majority of our 79 DZ have S/N ratios of 5-6. 
While this is sufficient to identify the pressure broadened ab¬ 
sorption features characteristic of cool DZ, detecting fields 
below 2 MG is not possible. For example, SDSS1152-I-1605 
was identified from a spectrum o bserved with the SDSS 
spectrogr aph ( Koester et al.| |201ll) at a S/N ratio of ~5. 
Although iKoester et al.l ( 201 il l speculated that this object 
may be magnetic, the quality of the data available to them 
was insufficient for a firm conlusion. Its newer BOSS spec¬ 
trum clearly reveals Zeeman splitting (S/N ratio of 9). We 
show a comparison of these two spectra in Figure [5] 

Inspecting the cumulative distribution in field strengths 
(Figure [S]) demonstrates the difficulty in detecting mag¬ 
netism for Bs < 2 MG. Above ~1.9 MG, the distribution is 
approximately linear in log (Bs), which i s consistent with th e 
distribution seen in other WD samples dKawka et al.|[2007l) . 
Below this value the only MWD found is SDSS 1330+3029 
(Bs = 0.65 MG), which is made possible by the exceptional 
S/N of its spectrum. This suggests that several objects with 
lower quality spectra may have magnetic fields between 0.65 
and 1.9 MG. 

As for resolution, the BOSS instrument has a resolving 


4 The dominant effects on the location in (u — g, g — r ) colour 
space are the effective temperature and the atmospheric abun¬ 
dance of Ca which causes strong opacity in the u-band. 



Wavelength [A] 


Figure 5. Mg triplets for both SDSS spectra of SDSS1152+1605. 
The Zeeman splitting on the earlier spectrum (left) is very am¬ 
biguous, yet well resolved in the newer BOSS spectrum (right). 



Figure 6. Cumulative total of magnetic objects from our sample 
versus the log of magnetic field strength. The increase is seen to 
be approximately linear in log(Bs) above 2 MG. 

power of approximately 2000 dSmee et al.ll2013l ). This im¬ 
plies a minimum detectable field via Zeeman splitting of a 
few hundred kG, and so SDSS1330+3029 is representative 
of the lowest detectable field for DZ in SDSS. 


5.2 Magnetic field origin and evolution 


While the true incidence of magnetism among WDs 
on the whole is still widely debated, estimates be¬ 
tween 5-10 percent are common for isola te d degenerate 
objec ts dWickramasin ghe fe Ferrarid l200d : iLiebert et aid 
120031 : ISion et ah 20141) . It has been suggested that older 
cooler) WDs exhibit a higher incidence of magnetism 
Kawka fe Vennesll2014l : [Liebert et al.ll2003l ). which at face 
value is supported by the large fraction of cool MWDs in 
our sample. However, we can at present not exclude that 
the high incidence of magnetism is linked to the presence 
of metals in the atmospheres of the cool DZ, e.g. through 
(merger) interaction with planets. 
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8 M.A. Hollands 


In addition, the origin of magnetic fields among WDs re¬ 
mains under discussion with the most plausible mechanisms 
proposed being: 


(i) From the intial-to-final mass relation (ICatalan et al.l 
120081') and main sequence lifetime as a function of stel¬ 
lar mass, most of the known WDs are thought to 
have evolved from A and B type s tars. These stars 
are known to exhibit magnetic Helds (lAngel et ahl Il98ll: 
IWickramasinghe fe Ferrario l2000l : iNeiner et al. 2014 ), with 
the peculiar Ap and Bp stars having comparatively higher 
fields. As the star evolves off of the main sequence it is 
expected that the magnetic flux of the progenitor star 
is conserved and so the change in stellar radius ampli¬ 
fies the surface field, i.e. Bwd/Bms = (-Rms/Rwd) 2 - 
This is known as the fossil field hypothesis, and can 
produce WDs with field strengths in the observed range 
(l\Voltierlll964l: lA ngel fe Landstre etlll970l : lAngel et al.lll98ll ; 
IWickramasinghe fe Ferrariol 2000lj . 

Ohmic decay is expected to cause magnetic fields to de¬ 
crease in strength with time. However, the timescale for 
this is expected to be of the order 10 10 yr due to the 
high electrical condu ctivity in the degenerate cores of WDs 
(IWendell et al.|[l987l 'l . Therefore the fossil field hypothesis is 
not unreasonable for describing the field origin in the old 
WDs we identify in this work. 

However, recent estimates of magnetic incidence ex¬ 
ceeding 10 perce nt dLiebert et al.l I 2 OO 3 I : iKawka et al.l 120071 : 
ISion et al.l 120141 1. challenge the fossil field hypothesis. The 
space density of Ap/Bp stars is insufficie nt to account for 
all the known MWDs with Bs > 1 MG dKawka fe Vennesl 
l2004f) . and so at least one other evolution channel is required 
for pr oducing MW Ds. 

(ii) I Tout et ahl (l2008l ~) suggested that WDs with Bs > 
1 MG are the products of an initial binary origin. Stellar 
evolution of one of the binary components can lead to a 
common envelope (CE) stage. It is during this phase that 
a magnetic dynamo may be generated within the CE. The 
resulting field then persists beyond the lifetime of the CE, 
within the now close binary or merged single object. For a 
close binary, a merger may take place later. 

The binary origin of these highly magnetised WDs nat¬ 
urally leads to higher masses than the canonical 0.6 Mq 
for non-magnetic WDs, compatible with the observation 
that MWDs are typically more m assive than non-magnetics 
dLiebertll 19881 : iLiebert et al.11200^ ') . 

However, a binary origin would in our case raise ques¬ 
tions about how these WDs come to be polluted by ma¬ 
terial from a remnant planetary system. This model need 
not be constrained only to stellar binaries. iNordhaus et aid 
d201l|) suggested that the engulfment of gaseous planets 
or brown dwarf companions during the asymptotic giant 
branch (AGB) phase could also lead to magnetic dynamo 
generation and eventually a high field mwd. in this sce¬ 
nario, mwds would not be expected to have higher masses 
than non-magnetics, but it would allow for evolved plan- 
etary systems wh ich later pollute the wd with metals. 
iFarihi et al.l (1201 il l identified a cool (T e s= 5310 k) magnetic 
(B ~ 120 kG) DAZ white dwarf, and speculate on that ba¬ 
sis that the WD underwent a CE with a closely orbiting gas 
giant planet during the progenitor star’s AGB phase, lead¬ 
ing to the emergence of a magnetic field. If this is indeed 


the mechanism from which magnetic fields are produced in 
DZ, it may explain the particularly high magnetic incidence 
found in our sample. 

Unlike the fossil field hypothesis, the giant planet CE sce¬ 
nario would be correlated with the presence of metals in 
the atmospheres of WDs, where the metal lines are an in¬ 
dicator of an evolved planetary system. Therefore, if DC 
white dwar fs, which originate from the same stellar popula¬ 
tion as DZ dFarihi et al.ll2010l j. have a significantly different 
distribution of magnetic fields, then this would present a 
compelling case for the CE hypothesis. 

(iii) An alternative origin for magnetism among WDs is 
aui dynamo generation. For a differentially rotating WD 
with a convective envelope, a magnetic dy namo may be 
gener ated at the base of the convection zone dMarkiel et al.l 
1 19941) . However this wou l d be unl ik ely to produce fields on 
the order of 1 MG dThomas et al.l Il995l ) , and would lead 
to magnetic fields strongly aligned wit h the WD ro tation 
axis which is in general not observed (lLatter et al.1 Il987l : 
iBurleigh et al.lll999l : lEuchner et al.ll20051 h 

5.3 The apparent lack of magnetism in warm DZ 

T he largest sample of WD s identified in SDSS was presented 
bv lKleinman et aD (120131 '), using SDSS DR7 spectroscopy. In 
total they identified 257 DZ, most of which are hotter than 
the sample we present here (T e ff> 8000 K), in which case Ca 
H/K are usually the only metal lines detected. Unlike in cool 
DZ where the broad wings of the H/K lines absorb most of 
the flux below 4000 A (Figure [5}, using these lines to detect 
> 1 MG fields should in principle be a trivial task. Addi¬ 
tionally, because of the larger sample size, there is also an 
abundance of these spectra where the H/K lines have good 
S/N ratios. We find 64 spectra with S/N > 10 (25percent), 
and 27 with S/N > 15 (10percent). 

Inspecting the Ca H/K lines of all 257 DZ did not 
reveal magnetic splitting for a single object. This is in 
stark contrast to our fraction of 13percent. One object, 
SDSSJ080131.15+532900.8, has what appear to be broad¬ 
ened Ca H/K lines which could indicate a magnetic field. 
However the SDSS images reveal this WD to be situated 
~7 arcsec away from a bright (r = 13.6) M star, which likely 
caused flux contamination in the DZ spectrum (obtained 
through a 3 arcsec fibre). 

The lack of any magnetic DZ in this sample either 
suggests a diffe r ent se t of selection effects at work in the 
iKleinman et al.l (|2013l i and our own, or that the incidence 
of magnetic fields in DZ differ s above and bel o w 90 00 K. 
The selection procedure used bv lKleinman et ail (l2013l l was 
distinctly different to the methodology we employed e m ■ 
They fitted DZ templates to a sample of WD candidate spec¬ 
tra, and so magnetic objects will have evaded detection if 
Zeeman split lines strongly affected the \ 2 °f their fits. 

We performed an independent check by inspect¬ 
ing the warm DZ identified from SDSS DR10 by 
ICentile Fusillo et al] (12014, submitted) . They selected can¬ 
didate WDs from a colour-cut in the (u — g ) vs. (g — r ) 
plane situated above the main sequence compared to our cut 
in Figure [l] and making use of proper-motions. All spectro¬ 
scopic objects with g < 19 bounded by this cut were visu¬ 
ally inspected and classified into the various WD subclasses. 
We inspected the H/K lines of all objects classed as DZ for 
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The incidence of magnetic fields in cool DZ white dwarfs 9 


splitting. Again, of the 118 unique objects, we did not find 
a single star t hat can be convincingly claimed to be mag¬ 
netic. Since the lGentile Fusillo et al.l (12014, submitted! ) sam¬ 
ple were only selected by colour and proper- motion, they are 
not bi ased against finding magnetic DZ as iKleinman et al.l 
J2013f) might be. 


Kepler et al.l (l2015h have recently published a list of new 


WDs from SDSS DR10 spectra including 397 objects clas¬ 
sified as DZ (where mos t have T e ff> 10000 K). Inspecting 
these 397 spectra reveals [Kepler et al.l d2015li have indepen¬ 
dently discovered 4 of the WDs that we have shown are mag¬ 
netic (SDSS0735+2057, SDSS0832+4109, SDSS1003-0031, 
and SDSS 1536+4205), however, beyond these, we found no 
further magnetic objects. 

The dearth of MWDs with T e g> 8000 K in the above 
3 samples may suggest that magnetic fields are preferen¬ 
tially generated several Gyr after leaving the main-sequence, 
at least for WDs with evolved planetary systems. If this is 
the case, then the origin of these Helds remains a difficult 
question to answer. We speculate that the gas giant/CE 
scenario (H5.2I1 previously proposed to o ccur on the AGB 
dNordhaus et al.ll201ll ; IFarihi et al.ll201lh . may still remain 
a viab le p ossibility at long W D cooling times. IVeras et al.l 
(l2013ll and lMustill et all (120141 ') have performed simulations 
indicating that planet-planet scattering can still occur many 
Gyr after stellar evolution to the WD stage, and show that 
these scattered planets will in some cases collide with cen¬ 
tral star. We suggest such a collision with a scattered gas 
giant might lead to magnetic field generation. However this 
scen ario has two prob lems : the small fraction of systems 
that IVeras et al.l (120131) and lMustill et all (l2014h expect this 
to occur for in comparison to our magnetic incidence lower 
limit of 13 ± 4 percent, and the lack of hydrogen lines ob¬ 
served in of the spectra in Figure [2] 


5.4 Comparison with magnetic DAZ 


For cool DZ we arrive at a magnetic incidence of 13 ± 
4 percent. Yet if we compare this against DAZ WDs, the 
result is very different. In fact very few magnetic DAZ are 
known at all, and their magnetic fields are not nearly as 
stro ng as found a m ong t he DZ in this study. 

Kawka et al.l (1200711 list all magnetic WDs known 
up to June 2006. Among these are the 3 previously 
known magnetic DZ (SDSS0157+0033, SDSS1214—0234, 
SDSS 1330+3029). However, not a single magnetic DAZ was 
known at that time. 

Since then, 4 m agnetic DAZ have been identified, 
all wi th T eff < 8000 K (IFarihi et alj 20 111 : IZuckerman et all 
l201ll : iKawka fe Vennesl 12011 . 2014), and with the most 
magnetic (NLTT 10480) possessing a field of only 0.5 MG 
dKawka fe Vennesll201li f . As with DZ ( 45.311 all known mag¬ 
netic DAZ have T e ff> 8000 K, again suggesting field gener¬ 
ation late on the WD cooling track. 

Additionally, the fact that cool DAZ are not found with 
the same regime of magnetic field strengths as DZ is some¬ 
what surprising as they will have similar cooling ages, and 
so the magnetic field distribution would be expected to be 
the same, assuming similar progenitors. 

One possible explanation is that because metal lines in 
DAZ appear weaker for a given metal abundance compared 
with DZ, magnetic spliting of these lines will smear them out 


in the continuum. Therefore strongly magnetic DAZ would 
instead be classified as magnetic DA (where the magnetic 
field can still be inferred from the Balmer series). Higher 
S/N spectra may reveal known magnetic DA WDs to also 
be metal polluted. 

It is also worth noting that, to date, there are no known 
magnetic DBZ. This is rather peculiar considering that met¬ 
als produce stronger lines in atmospheres dominated by he¬ 
lium than those dominated by hydrogen, which should easily 
allow the detection of fields up to a few MG. 


5.5 Field topology 

Our results from 1 )3.31 showed that an offset dipole topology 
provides a reasonable explanation for the minimal magnetic 
broadening seen in the Zeeman lines of some of the WDs 
shown in Figure [2] 

In this scenario, SDSS1536+4205 has a dipole offset 
away from the Earth leading to the distribution of observed 
field strengths appearing sharply peaked at 9.59 MG. It fol¬ 
lows from this model that the opposite, invisible hemisphere 
of the star exhibits a large gradient in field strengths with 
a strongly magnetic spot (~50MG) emerging at the pole. 
If the sign of a z was reversed, the strong gradient in the 
field across the visible hemisphere would have major obser¬ 
vational consequences. The er components in the Mg and 
Na triplets would be magnetically broadened to the extent 
of reducing their depth to only ~ 15 percent of the 7r compo¬ 
nent. Therefore, identifying the magnetic nature of the star 
would require a S/N ratio of at least 40 (for a 3<r detection). 

If the offset dipole model is the correct interpretation for 
the narrow Zeeman lines, this has a profound consequence 
for the incidence of magnetism in cool WDs. As discussed 
above, of these 10 MWDs in our sample, SDSS0157+0033, 
SDSS1003—0031, SDSS1536+4205, and SDSS2325+0448 all 
have Zeeman triplets that could arguably be explained by 
the offset dipole model. If cool MWDs have a tendency for 
their dipoles to be offset, and if all 4 of the above stars are 
viewed with their dipoles offset away from us, then statisti¬ 
cally this implies that several of the other WDs within our 
full sample should have dipoles offset towards the Earth. The 
o components of their Zeeman lines would be broadened to 
the point that they cannot be distinguished at the low S/N 
of the SDSS spectra, and so they would not been identi¬ 
fied as magnetic. Thus, if true, the offset dipole scenario 
increases the selection bias against identifying magnetism 
in cool WDs. 

An alternative explanation for the narrowness of the 
Zeeman lines could also come from a non-uniform distribu¬ 
tion of metals across the surface of the WD. If for instance 
the accreted material accumulated at the poles, the result¬ 
ing spectrum would exhibit splitting consistent with only 
the polar field strength. However, to reproduce the mini¬ 
mal magnetic broadening we observe, the metals would have 
to be constrained to such a small region that the resulting 
Zeeman triplets would have negligible depths. Alternatively, 
metals confined to the magnetic equator would also produce 
a spectrum showing a small range of field strengths, but 
would be able cover a much greater portion of the visible 
surface without significant magnetic broadening. However, 
equatorial accretion would necessitate at least a quadrapolar 
field. 


© 2015 RAS, MNRAS 000. HI fill 























































10 M.A. Hollands 


iMetzger et al.l d2012i 'l considered accretion of metals 
onto a MWD. If the sublimation radius of the WD is smaller 
than the Alfven radius (true for our MWD sample), then 
material is expected to accrete along the mag netic field lines 
as soo n as it enters an ionised state. However. IMetzger et al.l 
(l2012l 'l also discuss a potential caveat to this scenario. The 
presence of dust grains mixed within the gaseous disc may 
inhibit the ionisation of the gas component, and so even a 
strong magnetic field may have little influence over the ac¬ 
cretion flow of rocky debris. 

It is possible to distinguish between the two scenarios 
of non-uniform metal accretion and an of fset dipole by con¬ 
siderin g the rotation of magnetic WDs. iBrinkworth et all 
(l2013ll showed that isolated magnetic WDs have typical ro¬ 
tation periods of hours to days. Many of the WDs they ob¬ 
served in their sample have effective temperatures that fall 
into the range we study here. In general the ma gnetic axis of 
MWDs are not aligned with the rotation axis IILatter et al.l 
1 19871 : iBurleigh et all Il999l : lEuchner et al.l 120051 1 . therefore 
by taking spectra at multiple epochs, one would expect to 
see variation in the Zeeman line profiles of these stars. 

For non-uniformly distributed metals, the projected 
area of the metal-polluted region would phase with the rota¬ 
tion period. Spectroscopically, this would be seen as a reduc¬ 
tion in depth of all Zeeman lines in the triplet proportional 
to the change in projected area. 

For an offset dipole, the effect of rotation would be 
to bring more/less of the concealed, highly-magnetised pole 
into view, leading to increased/decreased broadening of the 
a components of the Zeeman triplet. However the n compo¬ 
nent would remain unchanged in depth. 


6 SUMMARY 

We have identified a sample of 79 DZ with T e g < 9000 K of 
which at least 10 possess magnetic fields in the range 0.5- 
10 MG. This implies a minimum incidence of 13 ± 4 percent 
which is substantially higher than for young hot DAs. Ac¬ 
counting for various sources of bias, such as poor signal-to- 
noise spectra and that these objects are identified from a 
magnitude limited sample, suggests that the true incidence 
of magnetism in cool DZ is almost certainly higher. We also 
note the narrowness of the Zeeman lines in several of these 
objects and show that these are most likely the result of 
a complex field topology such as an offset dipole. The si¬ 
multaneous release of SDSS DR11 and DR12 will provide a 
plethora of new DZ spectra from which improved statistics 
can be calculated. 
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